Mechanism for current saturation and energy dissipation in graphene transistors 
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From a combination of careful and detailed theoretical and experimental studies, we demonstrate 
that the Boltzmann theory including all scattering mechanisms gives an excellent account, with no 
adjustable parameters, of high electric field transport in single as well as double-oxide graphene 
transistors. We further show unambiguously that scattering from the substrate and superstrate 
surface optical (SO) phonons governs the high field transport and heat dissipation over a wide range 
of experimentally relevant parameters. Models that neglect SO phonons altogether or treat them 
in a simple phenomenological manner are inadequate. We outline possible strategies for achieving 
higher current and complete saturation in graphene devices. 



Electronic devices operating in the radio frequency 
(RF) regime play a central role in modern communica- 
tion. Graphene possesses a compelling potential in high- 
efficiency RF analog devices due to its excellent carrier 
mobility /i > 10,000 cm^/Vs at room temperature. Crit- 
ical to the operation of such devices is the realization and 
control of carrier velocity saturation at high source-drain 
electric field in graphene transistors. This subject has 
attracted numerous experimental and theoretical stud- 
ies recentl}^^"^, but no comprehensive analysis of experi- 
ments has been reported, which is the topic of this Letter. 

Graphene transistors operating in the small source- 
drain bias Vsd regime are described by the Drude model, 
where the current density j = neu = ne^E, depends 
on the electron density n, the transverse electric field 
E = Vsd/L and the mobility /i. At large electron drift 
velocity, inelastic collisions with phonons become increas- 
ingly frequent, eventually leading to a saturated Usat- It 
is in this velocity (current) saturated regime where ana- 
log amplifiers operate. Recent high-field transport mea- 
surements have demonstrated a high current density j of 
a few mA//xm in graphene transistors, with a few V on 
14d- This remarkable current-carrying capability is com- 
parable to that of carbon nanotubes (CNTs)2H2l and ex- 
ceeds the performance of silicon transistorgiSl. In contrast 
to CNTs, where the high-bias current in sufficiently long 
tubes reaches a full saturation of approximately 25 /xAPt^, 
the current in graphene does not yet fully saturate^El^ ex- 
cept in the presence of a carrier density gradienffl The 
saturation in CNTs can be described by a simple phe- 
nomenological model considering instantaneous emission 
of CNT optical phonons oi hui ^ 200 meV and zone- 
boundary phonons of hu « 160 meV^, although alter- 
native interpretations involving the SO phonons of the 
Si02 substrate are also plausible . For graphene, sev- 
eral phonon scattering mecha nisms, including the acous- 
tic phonons of gr aphen d^^ l ^^ l and the SO phonongi^ of 
the Si02 substrate^^'ISl have been shown to affect low- 
bias transport. In the high-bias regime, both optical 
phonons of graphene^El and the SO phonons of the Si02 
substrate^HI have been individually proposed and used 
to analyze experiments, although a comprehensive study 
combining all scattering channels and a careful evalua- 
tion of their individual contributions is still lacking. The 



closely related issue of dissipation of hot carrier energy 
becomes critical for devices operating in the high-curent 
regime. Recent modeling points to the crucial role of the 
Si02 substrate in graphene and CNT devices, although 
the detail s of the dissipation mechanism are yet to be 
articulatecPiSl. 

We obtain single-layer graphene sheet through me- 
chanical exfoliation and fabricate field effect transis- 
tors (FETs) on Si02 /doped Si substrates using standard 
lithographic techniques. Two-terminal and Hall bar con- 
figurations patterned onto a single graphene sheet allow 
us to correlate measurements of /sd(V^d)(TV), the small- 
bias resistivity p(V\^g) and its temperature dependence 
p{T) on the same sample (Fig.[l](a)). Measurements are 
carried out in a He^ cryostat at T=20 K unless other- 
wise noted. Fabrication and measurement details can be 
found in Ref s . [T7I and HM Samples used in this study show 
a low-density field effect mobility /ipE = {da/dn){l/e) of 
5,000 — 9,000 c m^/V s, comparable to samples used in 
other experimentgiEIll. 

A typical I-V trace is shown in Fig. IT] (e). Data 
used here are obtained from samples that showed negligi- 
ble change in p(Vbg) before and after I-V measurements 
(Fig. [l] (b)). I-V data are indepdent of the sweeping 
rate of Vgd, indicating that the sample reached an equili- 
brated state during the sweep. The two-terminal config- 
uration shown in Fig. [I] (a) faciliates a uniform current 
flow between the source and drain contacts and hence an 
accurate determination of j. It is crucial, however, to 
account for the contact resistance, -Rcon, in deducing the 
electric field E = {Vgd — IRcon)/L and the Joule heat 
P = I X (T4d — I Rcon)- Rcon is determined from mag- 
netotransport data, as shown in Fig. md), and lies in 
the range of 250-350 ft, where approximately 150 il is 
due to the cryostat wiring and several tens of fis may 
come from the patterned Au electrodes. Given its small- 
ness and ohmic nature at low biases, it is a reasonable 
assumption that the Au-graphene interface resistance re- 
mains constant during our I-V measurements. The drift 
velocity, u, is obtained directly from the current, /, us- 
ing u = j /ne and j = I /W where W is the width of our 
samples. 

Our theoretical treatment of high field transport in 
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FIG. 1. (Color Online) Field effect transistor device configu- 
rations and transport properties, (a) Optical micrograph of a 
device on Si02 /doped Si substrate. The graphene flake is out- 
lined in red. (b) Resistivity p{Vg) of a device before (black) 
and after (red) high-bias measurements, (c) Temperature- 
dependent resistivity p(T) at n = 1.30(2) x 10^^ cm"^. The 
coupling strength of LA (graphene) and SO (Si02 substrate) 
phonons are extracted from the fitting (solid line) . See Ref . [19] 
for details, (d) Two-terminal magnetoresistance 7?2pt vs mag- 
netic field. The contact resistance is -Rcon = 340f2, determined 
by extrapolating the magnetoresistance at the quantum hall 
plateaus, R2pt{B) = -Rcon -I- -§- (with the integer filling factor 
1/ shown on the figure), to S = 0. n= 1.82(2) x 10^^ cm~^ 
T = 20 K unless otherwise noted, (e) The measured (cir- 
cles) and calculated (lines) drift velocity (in units of Fermi 
velocity, vf) vs. electric field. The conduction of our devices 
tends to drop sharply and irreversibly near _E ~ 1 V/pm, 
presumably due to burning at local hot spots. Dashed line: 
theory with only impurities and the LA and LO phonons of 
graphene. Dash-dotted line: theory with only impurities and 
the SO phonons of the Si02 substrate. Solid line: the full 
theory including all scattering mechanisms. The density of 
holes is n = 2.09(2) x 10^^ 



3^^ cm~^, the density of charged im- 
purities is riimp ~ 5.8 X 10^^ cm~^. Experimental error bar 
is smaller than the size of the symbol. We fix the substrate 
temperature at rs=20 K and vary the graphene lattice temer- 
ature Tl with E from 20-380 K in this calculation (with Tl 
as high as 470 K for other measurements) . (f ) Comparison of 
experimental and theoretical drift velocity for several electron 
(right) and hole (left) densities a,t E — 0.6 V//im. 



graphene is based on the Boltzmann equation. 
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where the distribution function is assumed to be time 
independent and spatially uniform. The collision term 
on the right originates from impurities and phonons; 
We model the impurity contribution vifith charged and 
neut ral co mponents, labeled as imp and imp' respec- 
tively^SEIl. Phonons considered here are longitudinal 
acoustic (LA) and longitudinal optical (LO) phonons of 
graphen e, and t he surface optical (SO) phonons of the 
substrat6i^E3ESl_ The electron-electron interaction is in- 



cluded implicitly by choosing the "displaced" distribu- 
tion function, explained below. The collision integral for 
phonons, Sco\, has the form 
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f )<5(£fcQ — £p7 + shujq), with s — +1 {s — —1) is for 
phonon absorption (emission) and M^J is the matrix el- 
ement for the scattering process which takes an electron 
of momentum k in band a to momentum p in band 7. 
The quantities Nq and hojq denote the phonon occupa- 
tion factor and the phonon energy. We assume the equi- 
librium Bose-Einstein distribution for phonons, allowing 
for elevated temperatures for the graphene lattice and 
the substrate, but ignoring non-equilibrium hot phonon 
effects^S. The net charge density relative to the neutral 
Dirac point is given by 



nh 



= ffE[/kc-(l-/k.)] 



where 5 = 4 arises from the spin / valley degeneracy. 

It is estimated that for large electron densities, the 
electron-electron scattering time is sufficiently shorfPS 
that electrons come to equilibrium before any other scat- 
tering processes occur. The rapid establishment of an 
equilibrium electron distribution is approximated by as- 
suming a displaced Fermi-Dirac distribution for the elec- 
trons, 
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where a — ±1 denotes the conduction or valance band, 
Sfca = ctvpk is the energy spectrum of graphene {vp = 
10^ cm/s), /J,c is the chemical potential, Te is the tempera- 
ture of the electrons in the moving frame of reference, and 
u is the drift velocity. This approximation is justified a 
posteriori for the carrier densities reported here. Follow- 
ing Bistritzer and MacDonald^', we solve the equations 
for momentum loss rate, Q = — enE = — ^X^ka '^'^coi, 
and power dissipation, P — — enE • u — — gX^ka ^kaSco\: 
(which can be derived from the Boltzmann equation) for 
the three variables u, /io, Tc to obtain various transport 
coefhcients. 

The concentration of charged and neutral impurities, 
nimp and riimp' , are extracted from a global fit of the den- 
sity dependence of the low-tem perat ure resistivity p(Vbg) 
to the charged impurity modepi'^. Small discrepancies 
between the global fit and the actual data are corrected 
for by adjusting ni,„p slightly for each carrier density, 
resulting in nimp ^ 4.4-4.6 x 10^^ cm^-^ for holes and 
'T-imp ^ 5.3-5.8 X 10^^ cm~^ for electrons, the difference 
between the two attributable to the asymmetry of the 
conductivity. A few percent variation in Tijmp for elec- 
trons or for holes has negligible effect on the saturation 
velocity. The contributions of the neutral impurities to 
the residual resistivity are p'™P — 12.9 fl for holes and 
9.7 il for electrons. 



The temperature dependence of the low-bias resistiv- 
ity in fig. ff|c) is fit to a sum of three terms. The residual 
resistivity is due to impurities. The linear term is at- 
tributed to acoustic phonons33^ we find a deformation 
potential D — 18 eV, in agreement with other worlilS, 
We fit the non-linear T dependence of the resis tivity to 
the Bose-Einstein distribution for the SO phonons^^IiSIill, 
Parameters used to calculate the SO phonon frequen- 
cies are obtained from our measurements!^ and Ref. [3S1 
Using the linearized Boltzmann equation with the re- 
laxation time ansatz, and making the approximation 
that the electron-SO phonon scattering strength is mo- 
mentum independent, we find p^'-' = X]i(^/4e^)ff|(l + 
huJi/eF)N{uJi) in the low temperature, zero-bias limit, 
which we fit to the nonlinear part of the low-bias data to 
extract the coupling constants gi (i = 1,2 refers to the 
two SO phonon modes). The low-bias data show negli- 
gible dependence on the LO phonons, so we are unable 
to deduce this coupling parameter from experiment, and 
we use instead the theoretical coupling constants from 
Refs. [26|and|2Il As shown in Fig. [5] the dominance of 
SO phonons renders our results insensitive to the cou- 
pling strength to the graphene LO phonons. 

The results obtained from a numerical solution of the 
above equations are plotted in Fig. [2] for two representa- 
tive sets of lattice and substrate temperatures (Tl and 
Tg, respectively). Here we have taken Tl — Tg for sim- 
plicity to illustrate the key features of our theory. A 
more sophisticated treatment of Tl is implemented in 
Fig. [1] to compare to data and the effect of Tg is ex- 
amined in Ref. UHl The most striking result is that SO 
phonons are the principal scattering mechanism for high 
field transport: the drift velocity and the electron tem- 
perature obtained from the theory including all phonons 
are extremely well approximated in a model that retains 
only the SO phonons (and impurities), and more than 
~95% of the power dissipation also occurs directly into 
the SO phonons. A surprising feature is that the inclu- 
sion of SO phonons leads to an increase of the electron 
drift velocity, counter to the intuition that additional 
scattering mechanism should decrease it. The origin of 
this behavior lies in the fact that SO phonons also pro- 
vide an efficient route for energy dissipation, leading to 
a drastic drop of electron temperatures, which translates 
into higher saturation velocities. 

A detailed comparison between theory and experiment 
requires us to carefully consider the dependence of Tl on 
the applied field, E. We use an empirical relation of lat- 
tice temperature vs power dissipation, obtained from Ra- 
man spectroscopy recentljBl to estimate Tl as a function 
of E in our devices. Tl rangeaiSfrom 20 to 380 K for the 
device shown in Fig. fl] (e) . The accuracy of Tl , as well the 
coupling strength of the LO phonons in graphene, plays 
a minor role in our calculations due to the dominance of 
the SO phonons. In Fig. 1, we keep the substrate tem- 
perature at the bath temperature Tg = Tbath = 20 K as 
our calculation of the heat flow through the Si02 sub- 
strate indicateJiS an upper bound of Tg = 250 K, which 
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FIG. 2. (Color Online) The electron drift velocity, elec- 
tron temperature and power dissipation in graphene field 
effect transistors. Panels (a), (b), and (c) consider a low- 
temperature device, with Tl = Ts = 100 K, whereas the 
remaining panels assume Tl = Tg = 600 K, which corre- 
sponds to a typical room-temperature device with current 
induced heating. In panels (a), (b), (d), and (e), results 
are shown including all phonons (solid lines); LA and LO 
phonons (dash-dotted lines); LA phonons (dashed lines); and 
SO phonons (open circles). Panels (c) and (f) show total 
power dissipation (dashed line) and also the power dissipa- 
tion into the SO phonons (solid line). For all plots, we take 

density of charged impu- 



carrier density n = 2 x lO^'^ cm ^ 



nties rii, 



5 X 10 cm , and neglect neutral impurities. 



produces only a small deviation (< 4%) from results ob- 
tained with Tg = 20 K. There are no adjustable parame- 
ters in our calculations. 

Fig. fl] (e) shows that our theoretical results including 
either all phonons or only SO phonons agree with the 
measured u(E) for a sample of density n — 2.09 x 10^^ 
cm~^ to better than 4%. The agreement is somewhat 
worse for samples of lower densities, with ^15% dis- 
agreement at the smallest density shown in Fig. IT] (f). 
We attribute the worse agreement for small densities to 
complications arising from the formation of electron-hole 
puddles in real samples. We stress that a neglect of 
SO phonons results in substantial qualitative disagree- 
ment between theory and experiment. A simple estima- 
tion of u can be obtained by assuming instantaneous 
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FIG. 3. (Color Online) Saturation current density jsat of 
graphene on Si02. The saturation current density jsat vs. 
carrier density n for charged impurity densities, nimp = 
5 X 10" cm"^ (solid) and nimp = 5 x 10^° cm~^ (dashed). 
Neutral impurities are neglected, jaat is approximately linear 
for n > 5 X lO^'^ cm~^ and is slightly lower for the higher riimp. 
jsat is defined as j at i? = 2 V//im; as seen in the inset, all 
currents have reached saturation at this field. Inset: current 
vs. electric field j{E) for three representative carrier densi- 
ties, n = 5 X 10" (cyan), 2 x 10^^ (blue), and 1 x lO" cm"^ 
(green). Solid and dashed curves correspond to the higher 
and lower nimp respectively. All data in this figure are calcu- 
lated at the representative substrate and lattice temperatures 
Ts = 20 K and Tl = 500 K. 



emission of the relevant optical phonon,^* which has 
been demonstrated to be quite reasonable for high field 
transport in CNTs^. Such a treatment for graphene^ 
leads to Usat/vp — hujso/E-p, and results in a calculated 
I = V/{Rimp + V/Isat) (dotted line), which is -20% lower 
than the experimental data at E — 0.6 V//ini. 

Figure [3] shows a theoretical prediction of saturated 



current density vs. carrier density in graphene transistors 
fabricated on Si02 for two charged impurity densities. 
Both curves display a linear regime for rt > 5 x 10^^ cm^^ , 
where j reaches a few mA//^m. These predictions point 
to the prospect of high-performance graphene linear am- 
plifiers. As the sample quality improves, the current satu- 
ration will occur at lower electric field, allowing for exper- 
imental access and a greater operational range for these 
devices (inset). 

Finally, given the crucial role played by the SO 
phonons, it is natural to wonder if the saturation cur- 
rent may be enhanced by using another substrate^^'^Il 
or a double-oxide graphene transistor. We have studied 
theoretically Hf02, AI2O3 and Zr02 substrates in vac- 
uum/graphene/oxide structures, and find that, in spite 
of the large range of uj^o involved, the variation of the 
saturation velocity is less than 25%. We have also stud- 
ied in detail, both experimentally and theor etical ly, the 
double-oxide Hf02/graphene/Si02 structure^^ESl, xhe 
measured drift velocity is in excellent (4%) agreement 
with theory, and is again dominated by the SO phonons, 
but, surprisingly, is lower than the drift velocity for SO2 
substrate alonel^. 

In summary, by combining careful experimental and 
theoretical studies, we demonstrate that at high electric 
field, hot electrons in graphene lose energy predominately 
by emitting surface optical phonons of the substrate. The 
resulting current saturation can be accurately explained 
by a Boltzmann theory using experimentally obtained 
inputs. Cleaner samples are necessary to achieve full 
velocity saturation at accessible bias field strengths of 
w 1 Y/fim. 
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